FORMULAS FOR LYAPUNOV EXPONENTS 



ALEXANDRE T. BARAVIERA AND PEDRO DUARTE 

Abstract. We derive a series summation formula for the average logarithm norm of the 
action of a matrix on the projective space. This formula is shown to be useful to evaluate 
some Lyapunov exponents of random SL-matrix cocycles, which include a special class for 
which H. Furstenberg had provided an explicit integral formula. 

<N 

o 

1. Introduction 

Lyapunov exponents measure the rate of separation of nearby orbits in a given dynamical 
system. Linear cocycles form a class of dynamical systems where the study of Lyapunov 



exponents is an important and active subject. Roughly, a linear cocycle is a dynamical system 
on a vector bundle which acts linearly on fibers, over some fixed dynamics on the base of the 
vector bundle. In this article we shall deal with a special class of linear cocycles where the 
base dynamics is a Bernoulli shift equipped with a (constant factor) product measure, and 
the fiber is the special linear group SL(d,K). These cocycles, that we will refer as 'random 
linear cocycles', were extensively studied by H. Furstenberg. See [8], [9], [10], [llj . More 
precisely, if G is some matrix Lie group, like SL(d, R), a random linear cocycle is determined 
by a probability measure [i on G together with a map F : G N x M d — > G N x M d defined by 
J> F(g, v) = {cr(g),go v), where g denotes a matrix sequence {g n }neN and a : G N —¥ G N denotes 

<^ the shift map o-{g n }n£N = {<7n+i}neN- The measure /ionG determines the product Bernoulli 

measure \i = /U N on the bundle's base G N , which is shift invariant. It is always assumed in this 
theory that \x is 'integrable', which means that J^log -1 " \\g\\ dfi(g) < +oo. The iterates of F 
are given explicitly by F n (g,v) = (a n g,g n -i . . . gigov). From a probabilistic point of view, 
this random cocycle is determined by the independent and identically distributed random 
process X n : G N — > G, defined by X n (g) = g n = Xo(a n g). Furstenberg and Kesten [10] have 
proven that the following limit, expressing the growth rate of the product random process 
^ Y n = X n _i . . . X\ Xq, always exists and is constant for //-almost every g £ G N , 

"X 

?— i 1 1 n 

d A(/u) = lim - log||cfa_i ... gig \\ = hm - log F n (#) . 

n— >+oo n n— >+oo n ~ 
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The average growth rate A(//) relates to the Lyapunov exponent as follows. Given (g, v) G 
qH x ]g>c^ Ly a p Unov exponent at g along v is the limit 

\^(v) = lim - log \\g n -i ■■■ 9i9o v\\ = lim - log \\Y n (g) v\\ , 

n— ¥+00 n n— Y+oo n " — " 



which exists for every v & M. — {0} and //-almost every g G G . Moreover, this limit is 
constant //-almost everywhere and \^(v) < A(//) for every v E M rf — {0}, with equality for 
almost every vGlf In fact, strict inequality can only occur for vectors in some //-invariant 
proper vector subspace V C M. d , i.e., one which is invariant under all matrices g in the 
support of //. See theorem 3.5 of [11]. This shows that A(/f) is the largest Lyapunov exponent 
of the random cocycle. Furstenberg and Kifer give a nice variational characterization of all 
the Lyapunov spectra in [11] . but we shall only deal with the largest Lyapunov exponent 
A(/i) here. In [9] H. Furstenberg established the following integral formula for the Lyapunov 
exponent 



= / / l°g||0z|| du(x)dn(g) , (1.1) 

JSL(d,R) Jpd- 1 

where P d_1 denotes the projective space of lines in M. d , and v stands for any maximal /it- 
stationary measure. A measure v on is said to be ^-stationary iff /i x v is an F- 
invariant measure. This amounts to say that v is a fixed point of the convolution operator 
Pfi,(v) = H * v = f 9*vdfj,(g), induced by /i on the space of Borel probability measures on 
P d_1 . A //-stationary measure is said to be maximal if it maximizes the left-hand-side integral 



in (1.1). See for instance section 3 of [llj for proofs of these facts. Furstenberg also found 
very general sufficient conditions for the largest Lyapunov exponent to be strictly positive. It 
is enough that the group G generated by the support of /i is non compact, and no subgroup of 
G with finite index is reducible. A group G generated by matrices in the support of \i is said 
to be reducible iff there is a non trivial decomposition of M. d as a direct sum of //-invariant 
subspaces of M. d . See theorem 8.6 in [9]. 



Furstenberg's formula (1.1) indicates a way of computing Lyapunov exponents. But still 



a couple of problems persists. 

(1) To determine the //-stationary measures explicitly. 

(2) To compute the following integral numerically 

R u (g) = [ \og\\gx\\ dv{x) . (1.2) 

In [9 J Furstenberg solves the first problem for a special class of measures on SL(d, M). See 
theorem 7.3 of [S]. In this paper we address the second problem mainly. First we consider 
the uniform Riemannian probability measure m on IP^ -1 and prove in section 3 that 
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Theorem A Given g G SL(d,K) with singular values < Ai < A2 < . . . Xd, if X* > Xd/V% 
then the following series converges absolutely 

00 1 / \ 2 \ ri / \ 2 \ T ' d 

r=l nH hr d =r ^ * ' ^ * ' 

where 



r! (2n - l)!!---(2r d - 1)!! 
n!--T d ! d(d + 2)---(d + 2r-2) 

Moreover, the coefficients ®rt]..-,r d form a permutation invariant probability distribution on 
the finite set 3 r = { (n, . . . , r^) G N d : r% + . . . + = r }. 



In the last section 4, we provide some applications of Theorem A. First, we give an explicit 
formula for the largest Lyapunov exponent of the random cocycles where Furstenberg was 
able to give explicit stationary measures. This formula is given in terms of the integrals 
Rm{g)i t° which we can apply Theorem A above. In a few special cases, these formulas 
are used in numerical computations of some Lyapunov exponents. A second motivation for 
proving Theorem A was the role played by the integral R m (g) in the following conjecture. 
For any dimension d > 2 and every g G SL(d,R), 

( log p{kg)dk>R m {g), (1.4) 

JSO{d) 

where p stands for the spectral radius and dk represents integration with respect to the 
normalized Haar measure in the special orthogonal group SO(d). This is conjectured in J6j 
Question 6.6]. An analogous result is proved in [7] for the unitary group in GL(d, C). Theorem 
A is based on the following more general result, to be proved in section 2. 

Theorem B Given a probability measure v G J'(P rf_1 ), and g G SL(d, M.) with singular values 
< Ai < A2 < • • • Xd, if A* > Xd/V2 then the following series converges absolutely 



^) = io g A,-E^ E •••(1-3) 

P=l Zr r 1+ -+r d =r V A * / V A * ' 



where 

eg..,r d (fc^)= , / Qn,..,r d (A:x)dKx), 

~i- " ^d- Jpd-i 

Qn,...,r d '■ P d_1 — > R stands for the function Q ri ,...,r d (xi, ■ ■ ■ ,Xd) = a; i n • • • ^i^? ari ^ k G 
0(d, R) is an orthogonal matrix such that kg T gk~ l is a diagonal. 

Moreover, the coefficients @rt]...,r d {k,v) form a probability distribution on the finite set 3 r = 
{(?!,..., r d ) GN d : r x + . . . + r d = r }. 

Let us remark that in theorem A the coefficients ®rfj...,r d ' 
(1) are given explicitly, 
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(2) do not depend on the singular value decomposition of g, and 

(3) are invariant under permutations of the indices (r±, . . . , r^). 
Although more general, all these properties fail in theorem B. 

Defining &ri,...,r d (v) := Qrt]...,r d {id,v), we have &rfj...,r d (k,i') = ®ru...,r d (,k u) , and the 
measure v = kv is /I-stationary with Jl = k/sk^ 1 . Hence, there is no loss of generality in 
assuming that k = id, i.e., g is a diagonal matrix, in theorem B. 

If we have upper bounds on some of the coefficients Qri}...,r d (kv) then we can use theorem 
B to obtain a corresponding lower bound for R v (g). For instance, in proposition [2] we shall 
derive the following universal upper bound 

V ©, ri (k, v) ( i - ^Pf . . . ( i - < ( i - f 



riH Vr d =r v *' v *' v * 

Then, combining this with theorem B, we can derive the trivial lower bound Rv{g) > logAi, 
where Ai denotes the least singular value of g. 

The family of functions Q ri rd : P d_1 — y M, with (ri, . . . , r^) £ N d , separates points in 
P d_1 . Hence, by Stone- Weirestrass' theorem, the linear space spanned by these monomials is 
a dense subalgebra of C(P d_1 ). In particular, the measure v is completely determined by the 
'momenta' ©n,... rAy)- If we could devise some convergent iterative scheme to approximate 
these '/U-stationary momenta', instead of the fx-stationary measure v, then we would apply 
theorem B and get bounds on the Lyapounov exponent R u (g). 



2. A General Formula 
Given integers r\ > 0, . . . , ra > 0, consider the function Q ri ,...,r d '■ P d_1 — > R 

Qr\,...,r d {x\, • • • , = X-y 1 . . . d . (2-1) 

This is a bounded function taking values between and 1. Setting r = r% + . . . + r<f, the 
minimum value of Q ri ,...,r d is while the maximum value, (ri/r) n ... (rd/r) r ' d , is attained 

at the projective points with coordinates ^i-y/Vi/r, . . . , ±y/rd/rj . Let J'(P d_1 ) denote the 

space of Borel probability measures on P d_1 . Throughout this section, v will denote any 
probability measure in y(P rf_1 ). 

Proposition 1. Let A £ gl(d, R) 6e a symmetric matrix of the form A = k~ 1 Dk, where 
k -1 is an orthogonal matrix consisting of A's eigenvectors and D = diag(Ai, . . . , Ad) is the 
corresponding eigenvalue matrix. Then 

(Ax,x) r du(x) = © S..,^) A ? • • • A ? , (2-2) 

riH \-r d =r 
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where 

eg..,r d (fc^)= , / Q n ,...,r d (kx)du(x). (2.3) 

Proof. Using the multinomial formula we get 

r 



ri \r d 
d ' 



f li>2^iXi) du(xi,...,x d )= V — r-^ — :(/ Q ri ,...,r q du) XY ... A 

which is the stated formula with A = D = diag(Ai, . . . , A<j). The general case follows in the 
same way but replacing x by kx before integrating with respect to v. Note that (Ax,x) = 
Y^i=i (ki x ) 2 with k x = (ki x, . . . , kd x) . □ 



Given a probability measure v G CP(P ), and integers n > 0, . . . , r d > we shall denote 



the matrix function defined in (2.3) by GJ^ r : 0(<i, K) — )■ M. 
Corollary 1. {B^ r } is a family of non-negative bounded functions such that 

@ri,...,r d (^) = 1 / or every matrix k £ 0(d,M) . 



riH hr d =r 

Proof. First 



o < e^..., rd < - (n/r)^ • • • (r d /ry« = -^...^L 

L ' ' d r\\ . . . r d l r r r\\ r^! 

because < Q ri ,...,r d < (ri/r) ri . . . (rd/r) Td . By proposition [TJ 

Yl n,...,r d W= / . , —} ^ —) Q ri ,...,r d (kx)du(x) 

n+-+r d =r ri +-+r d =r J¥ l ' ' " d ' 



riH hr-d=r 

(cc, x) r dk*v(x) = 1 . 

1 

□ 



Proof of theorem B. We use the following Taylor' series for the logarithm function 

°° 1 

logx = logrro - ^ - (l - Xq 1 x) r . 

r=l T 

Take xo = X 2 and b = I — \* 2 g T g. Note that Ar^ ||<?£|| 2 = {X* 2 g T gx,x) and 
b = k^ 1 diag (l — p-, . . . , 1 — k. Hence, by proposition 1 
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/ log \\gx\\ dv{x) = - I \og\\g x\\ 2 dv{x) 

Jpd-1 2 Jpd-1 



logA*-V— / {\-\- 2 \\gxf) r dv{x) 
^ 2r Jpd-i V J 

00 1 f 
logA*-V] — / (bx,x) r du(x) 

r=l riH Yr d =r ^ * ' 



2\ r l / \2\ r d 



riH \-r d 

The assumption A* > \d/V% implies that 

A? 



a = max 

Ki<d 



< 1 . 



The r.h.s. converges absolutely because the absolute value series is majorated by 

00 1 00 r 

E^ E @S.,.(mk = E 

r=l riH hr'd=r' r=l 



□ 



Remark 1. Assuming A* > y (Af + A^)/2 we /zaue 

i2\ \2 



a = max 



6(0,1). 



We end this section with one more remark. 



Proposition 2. Given g £ SL(d, M) with singular values < Ai < ... < A^ ; let k be an 
orthogonal matrix such that g T g = k^ 1 diag(Af, . . . , A^) k. If A* > <J {\\ + A^) /2 then 



riH hr d =r v * 7 v * 7 v * ' 

Proof. Combine corollary [T] with remark [T] □ 
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3. Spherical Integrals 
Theorem A follows from theorem B using next formula. 
Proposition 3. Given integers r\ > 0, . . . , r& > 0, 

(2r 1 -l)!!...(2r ri -l)!! 
Qn rddm ~ d{d + 2) ... (d + 2r-2) ' 1 j 

This formula involves the concept of double factorial, which relates with Euler's Gamma 
function. The double factorial is the recursive function defined over the natural numbers by 
the relation n\\ = n{n — 2)!! with initial conditions 0!! = (—1)!! = 1. The Gamma function, 
defined by the improper integral 

T(x) = \ t x - x e- l dt (x > 0) , 
Jo 

is a solution of the functional equation 

r(x + l) =xT(x) . (3.2) 

Since T(l) = e~ t dt = 1 it follows at once that Tin) = (n — 1)! for every n G N. In other 
words, the Gamma function is a real analytic interpolation of the usual factorial function 
over the natural numbers. Likewise, because T(l/2) = y/2 it follows easily by induction that 
for every n 6 N, 

1\ (2ra-l)!! r- , . 

n+- 2 ) = { —^^. (3.3) 

We refer pQ for a comprehensive treatment on the Gamma function. See formula (1.1.22) 
there for a justification of the value T(l/2) = y/2. 

The Gamma function can be used to provide explicit formulas for the volumes of spheres 
and balls. Let D d = {xGlR a! : ||x|| 2 <l}be the Euclidean unit disk and denote its volume 
by Vd- As above, let S d_1 be the Euclidean unit sphere, i.e., the boundary of D rf , and denote 
its area by A^_! = Ld-i Ida, where a stands for the measure induced by the canonical 
Euclidean induced metric on S d_1 . The Divergence theorem, together with a simple change 
of variables, may be used to establish the following recursive relations between these volumes 
(see appendix A of [1]) 

V d = ^V d _ 2 and A d _ 1 = dV d . (3.4) 
From these relations we deduce explicit formulas for the volumes of balls and spheres: 

= —. —-T and Aw_i = 2 — ,, . . (3.5) 

r(i + d/2) r(d/2) v ' 



To see this set = ir d l 2 /T(l + d/2). The functional equation (3.2) implies that satisfies 
the same recursive equation as V<2, 

7T d / 2 7T ^ rf - 2 )/ 2 _27T 

d ~ r(i + d/2) ~ d/2 r(i + (d - 2)/2) ~ d d ~ 2 ' 
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But since T(l/2) = y/2, it follows that Ui = 2 = Vi. Also Uo = 1 = Vo- Hence the equality 
V d = U d holds for all d > 0. Finally, by flO) we get 



A d _i = d\ d = d 



7T 



d/2 



7T 



ef/2 



r(l + d/2) r(d/2) 



Proof of proposition [3| We are going to reduce the integrals (3.1) to the following 
family of integrals introduced in [5] 



,rd) 



_2ri 



x^ rd (ixi . . . dx d 



Using the Divergence Theorem the author deduces a recurrence formula from which he gets 
the following explicit formula for every d > 2, and every r± > 0, . . . , > 0, 



id(n, ■■■,r d ) 



T(r 1 + \) ...T(r d +\) 

r(n + --- + r d + i + f; 



(3.6) 



Check formula (8) of [5J. It is also easy to see with a change of variables' argument that for 
any r-homogeneous function / : W 1 — > R (see Corollary 1 of [5]), 

1 



, x d ) dx\ 



dx d 



d + r 



f(xx,...,x d )da(xi,...,Xd). (3.7) 



Now, combining (3.7), (3.6), (3.5) and (3.3) we get 

^ n ...^ rd C?C7(xi,...,X rf ) 



1 

d + 2r 



Id(ri, ■ • - ,r d ) 



A d _i 

<i + 2r r(ri + i) ...r(r d + i) 



2^ 



d/2 



r(r + i + f) 



rp) *v ^ 1 2^ 
d + 2r (2n - 1)!! ... {2r d - 1)!! 
2 r+1 (r + f)(r-l + |) ... | 
(2n -l)!!... (2r d -l)!! 
d(d + 2) ... (d + 2r-2) 



□ 



FORMULAS FOR LYAPUNOV EXPONENTS 



9 



Proof of theorem A. In view of theorem B we just have to compute: 
>(<*) 



Q%l..,r d ( k i m ) = ri , r " , / , . Qr u ...,r d (kx)dm(x) 



<r* I P <r* I 



Qr 1 ,...,r d (x)dhm(x) = — r / Q ri ,...,r d dm 



ri! . . .r d ! y P d-i ri! . . .r d \ 

the last equality because k*m = m, for every orthogonal matrix k E 0(d,R). Combining 
this computation with proposition [3] we obtain formula (1.3) in theorem A. The coefficients 
eg... t r d are obviously positive rational numbers, whic h by corollary hj form a probability 
distribution on the set 3 r . An inspection to formula ( 1.3 ) shows these coefficients are invariant 
under permutations, i.e., &ri,...,r d = &r^ 1 ,...,r 7Td , for every permutation ir of {1, . . . , d}. □ 

4. Some Applications 

Denote by M m the space of probability measures (i in SL(d, R) that have m as ^-stationary 
measure, i.e., \i * m = m. The class M m is closed under orthogonal averages, i.e., if \i E M m 
then J SO (rfM) Kfidm(k) E M. 

Proposition 4. For any measure ^ E M m , ite Lyapunov exponent is 



A(/x) = / -Rm(s) d/4s0 

JSL(d,M) 



Proof. Follows from Furstenberg integral formula ( 1.1 ). □ 



Theorem A can then be used to approximate this Lyapunov exponent. A class of examples 
in M m are the so called orthogonally invariant measures. A probability [x E CP(SL(d, R)) is 
said to be orthogonally invariant if k^/i = fj, for every orthogonal matrix k E SO(d, R). We 
list some equivalent characterizations of orthogonally invariant measures. 

Proposition 5. Given a measure [i E IP(SL(d, R)), the following are equivalent: 

(1) [i is orthogonally invariant, 

(2) fi*5 p = m, \/p£F d - 1 , 

(3) fi* v = m, W E ^(P^ 1 ), 

(4) fi = rriK * 6, for some measure 6 E T(SL(o!, R)), 

where mx stands for the normalized Haar measure on K = SO(d, R). 

Proof. The proof is straightforward. □ 
Given a matrix g E SL(d, R), consider the measure 



^ = m K *S g = / 5 kg dm K (k) . (4.1) 

JsO(d,M.) 
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Proposition 6. The measure (4-1) is orthogonally invariant, and its Lyapunov exponent is 
A(/i) = R m {g). 

Proof. Since fi is orthogonally invariant we have fi* m = m, and hence by proposition [4] 



X(n) = / Rm(g') dfj,(') = / R m (kg)dm K (k) = R m (g) , 
JSh JSO 

because all matrices k g have the same singular values. □ 

Consider now the matrix family 

gt=( ^ t -i Id ) G SL(2d,M) (t > 1) , (4.2) 

where Id denotes the identity d x d matrix. Next proposition refers to the following orthog- 
onally invariant measure fit = mx * o~ gt ■ 

Proposition 7. For every t > 1, the Lyapunov exponent of fit is 

~ 1 d(d + 2)...(d + 2r-2) ( 1 
XM ~ ^2r (2d)(2d + 2)...(2d + 2r-2) ^ " ¥ 

Proof. By proposition [6j \2d{^t) = Rm(gt)- Notice that matrix gt has d singular values 
equal to t > 1, and d singular values equal to 1/t < 1. Thus, applying theorem A with A* = t 



A M (A*t) = #m(<7t) = logt - £ ^ e ?. < ?,0,n,...,r <J f 1 " ^ 

r=l ri+. ,.+r [ j=r ^ 

Notice that if + . . . + + n + . . . + r d = r, 

R (2d) _ H (2n-l)!!---(2r fc -l)!! 

o,..,o,n,..,r d ri !... rd ! (2d)(2d + 2)---(2d + 2r-2) 



r 



(d) (2rf)(2rf + 2)---(2d + 2r-2) 
ri '-' rd d(d + 2)-..(d + 2r-2) 



Hence, because E n +...+r d =r & rf,...,r d = 1, 

R ( 2 «0 (2tQ(2rf + 2)---(2rf + 2r-2) 

2^ w o ) ...,o,n,..,r t( d(d + 2)---(d + 2r-2) ' 

ri+...+r ( ;=r \ / * ; 

and we get the given formula for X-zdivt)- n 

Corollary 2. T/ie function sequence A2d(A*t) increases with d, and for every t>\ 

1 / 1 + i 4 
lim X 2 d{vt) = logi + - log 



d ^ +00 zau ^ to 1 2 to V 2t 4 
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Proof. Notice that 

d(d + 2)...(d + 2r-2) 1 

(2d)(2d + 2)...(2d + 2r-2) ~ T ' 

and the right hand side decreases to 2~ r as d grows to +00. The series Y^T=i 2 r + 1 r ~~ f) 
converges absolutely and uniformly to the function 

^)= 1 °g*-E^v( 1 -?J • 

r=l v 7 

Because this series is essentially a geometric one we can compute its sum explicitly 

g(t) = logt + - log 

Then, by Lebesgue monotone convergence theorem lim^oo A2d(/^t) = <?(£)■ □ 




Figure 1. The graphs of t h-> \2d(pt), for the half-dimensions d = 1, 2, 3, 4. 

This corollary shows that for large dimensions, X^di^t) ~ l°g^ = l°g ||fft|l> which is somehow 
expectable since all matrices in the support of \xt have norm t. 

The graphs of these functions, computed in Mathematica are depicted in figure [TJ The 
dashed line represents the graph of g(t). 

For the following class of measures Furstenberg was able to give explicit stationary mea- 
sures, see theorem 7.3 of [9j. Given two probability measures \i\ and ^2 in SL(d,K), define 
the measure 

A* = Mi * m K * = / / / S gikg2 dni(gi)dm K (k)diJ,2{g2) ■ (4.3) 

JSL ./SO JSL 
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The //-stationary measure of is \x\ * m. In fact, by item 4. of proposition [5j the measure 
itik * M2 * Ml is orthogonally invariant. Hence (mx * M2 * Mi) * m = m an d 

/i * (/ii * m) = (/ii * rrix * M2) * (mi * Tn) 
= fii * (rriK * fJ-2* Mi) * Tn 
= * m , 

which shows that /ii * m is /U-stationary. 



Proposition 8. The Lyapunov exponent of (4-3) is 

A(m) 

Proof. By Furstenberg formula 



SL JSL 



a(m) 



SL JSL JSO 



SL JSL JSL JSO 



log 1 1 gife^a; 1 1 d(Mi * m)(sc) dm K (k) dfii(gi) 4*2(52) 

dm(x) 4*1 (tf'i) dm K (k) 4*1(01) 4*2(52) 



/ log 


5i«52,| / I, 


/pd-1 


5i x 



SL JSL JSL JSO 



Rrn(gikg2g'i) - R m (g[) dm K (k) 4*1 (ffi) 4*1 (51) 4*2(02) 
Rm(92g[) + Rg 2 g' m (gik) - R m {g[) dm K {k) 4*i(0i) 4^1(01) 4/2(02) 

SL ./SL ./SL JSO 

Rm(g2g[) - R m {g[) 4*1(01) 4*2(02) + 



SL JSL JSL V^SO 



SL JSL 

R g2 g' im {9lk) dm K (k) \ 4*i(0i) 4*i(0l) 4*2(02) 
#m(020i) 4*i(0i) 4*2(02) - / R m (g[) dni(g[) + / ^(01)4*1(01) 

SL ./SL JSL JSL 

Rm(929i) 4*i(0i) 4*2(02) ■ 



SL JSL 



On the fourth step we use item (b) of lemma [TJ and on the sixth step we use lemma [2j □ 



Lemma 1. Given g',g £ SL(d, M), 

(a) ^(5) < log 

(b) R m (g'g) = R m {g) + R gm (g'), 
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Proof. The proof of (a) is straightforward. Item (b) holds because 



Rgm(g') = / log \\g' x\\ dgm{x) = / log 



g 



\gx 

log \\g'gx\\ dm (x) — j log \\gx\\ dm(x) 

R m {g'g) - R m (g) 



Lemma 2. Given g G SL(d, M) and any probability measure v £ 'P(P d 1 ), 

R u (gk) dm K {k) = R m {g) . 



dm(x) 



□ 



/SO(d,K) 

Proof. The measure is orthogonally invariant. This because mx = mx * d~i, where / 
denotes the identity in SL(d, M), by item 4. of proposition [5j Then, by item 2. of the same 
proposition, mx * 5 X = m. Hence 

R u (gk) dmK{k) = / / log \\gkx\\ dis(x) dm,K(k) 
so Jso JP 1 *- 1 

/ log\\gkx\\ drnxik) du(x) 
P d - X Jso 



log ||<7z|| d{rriK * S x )(z) dm(x) 

l 

log \\gz\\ dm{z) dv(x) 

l 

R m (g) dv(x) = R m {g) ■ 



□ 



We consider now a special subclass of the previous. Given two matrices g',g € SL(d, R) 
define the measure 

/V ,9 = <V * m ^ * &9 = / <V% dm K (k) . (4.4) 

JSO(d,M) 

Corollary 3. TAie Lyapunov exponent of the measure fj, g ' !g is 

KVg^g) = Rm(g'g) ■ 

For example, A( / u g -i g ) = R m (I) = but the measure Ay-i, 9 is supported on a compact 
group g^ 1 SO(d, M) g, and hence should have zero Lyapunov exponent. 



Consider now the measure fi gs ,g t , where g s ,gt are matrices as defined in (4.2). By corollary 
HI KVg s ,gt) = Rm{g s gt) = Rm(g s t)- Hence 



1 , / 2s 4 i 



4 + 4 



Corollary 4. lira \ 2 d(Vg s ,g t ) = log(si) - - log . 

a— >+oo Z \ 1 + St 

Notice that log(si) = maxjlog \\gx\\ : g G supp(/x 5ajgt ) , x € P d_1 }, the norm of matrices 
in the support of fJ> ga , gt is not constant and A(// SsiSt ) is some kind of average of the logarithms 
log \\gx\\, with g 6 supp(^ 9sj9t ) and x G P rf_1 . Prom this we conclude that large dimensions 
bring the average \(fi ga . gt ) closer to its maximum possible value, log(si), provided st is large. 



A similar conclusion, assuming conjecture (1.4) to hold, is that for large t > 1 and large 
dimension d, 

log t = log \\g t || > / log p(kg t ) dm K (k) > R m (g t ) ~ log* . 
•/ SO 

Again, this shows that large dimensions bring the average f so log p(k gt) drnxik) close to the 
maximum value log ||<7t||. 



Acknowledgements 

This work was partially supported by Fundagao para a Ciencia e a Tecnologia through the 
project "Randomness in Deterministic Dynamical Systems and Applications" ref. 
PTDC/MAT/105448/2008. 



References 

[1] G. Andrews, R. Askey, R. Roy, Special Functions. Encycolpedia of Mathematics and its Appli- 
cations, (1999) Cambridge University Press. 

[2] A. Avila, J. Bochi, Lyapunov Exponents. Lecture Notes, School and Workshop on Dynamical 
Systems, ICTP-Trieste 

[3] A. Avila, J. Bochi, A formula with some applications to the theory of Lyapunov exponents. 

Israel J. Math., V. 131, (2002), pp 125-137 
[4] S. Axler, P. Bourdon, W. Ramey, Harmonic Function Theory , Graduate Texts in Mathematics 

137, (2001) Springer Verlag 
[5] J. Baker Integration Over Spheres and the Divergence Theorem for Balls. The American Math- 
ematical Monthly, Vol. 104, No. 1. (Jan., 1997), pp. 36-47. 
[6] K. Burns, C. Pugh, M. Shub and A. Wilkinson, Recents results about stable ergodicity, J. Stat. 

Phys. 113 (2003), pp 85-149. 
[7] J. P. Dedieu and M. Shub, On random and mean exponents for unitarily invariant probability 

measures on GL(n,C), in "Geometric Methods in Dynamical Systems (Il)-Volume in honor of 

Jacob Palis", Asterisque 287 (2003), pp 1-18. 
[8] H. Furstenberg, A Poisson Formula for Semi-Simple Lie Groups. Annales of Mathematics, Second 

Series, Vol. 77, No. 2, (1963), pp 335-386 
[9] H. Furstenberg, Noncommuting Random Products. Transactions of the American Mathematical 

Society, Vol. 108, No. 3, (1963), pp 377-428 
[10] H. Furstenberg, H. Kesten, Products of Random Matrices, The Annals of Mathematical Statistics 

Vol. 31, No. 2 (1960), pp. 457-469. 
[11] H. Furstenberg and Y. Kifer, Random Matrix Products And Measures On Projective Spaces, Israel 

Journal Of Mathematics Vol. 46, (1983) pp. 19-20. 



FORMULAS FOR LYAPUNOV EXPONENTS 



15 



Departamento de Matematica, Universidade Federal do Rio Grande do Sul, Porto Alegre, 
RS, Brasil 

E-mail address: baravi@mat.ufrgs.br 

Departamento de Matematica and Cmaf, Faculdade de Ciencias, Universidade de Lisboa, 
Campo Grande, Edificio C6, Piso 2, 1749-016 Lisboa, Portugal 
E-mail address: pduarte@ptmat.fc.ul.pt 



